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Abstract. We present the Ba even-to-odd isotope abundance ratios in 25 cool dwarf stars with the metallicity [Fe/H] ranging 
between 0.25 and -1.35. Our method takes advantage of the hyperfine structure (HFS) affecting the Ban resonance line of the 
odd isotopes. The fractional abundance of the odd isotopes of Ba is derived from a requirement that Ba abundances from the 
resonance line A4554 and subordinate lines /1 5853 and /1 6496 must be equal. The results are based on NLTE line formation 
and analysis of high resolution (R ~ 60000) high signal-to-noise (S/N > 200) observed spectra. We find that the fraction of the 
odd isotopes of Ba grows toward the lower Ba abundance (or metallicity) and the mean value in the thick disk stars equals 33 
± 4%. This indicates the higher contribution of the r-process to barium in the thick disk stars compared to the solar system 
matter The obtained fraction increases with the [Eu/Ba] abundance ratio growth in agreement with expectations. A significant 
fraction of the even isotopes of Ba found in old Galactic stars (the thick disk stars), ~ 67%, is in contrast to the prediction of 
the "classical" model of the .s-process and favors the value predicted by the "stellar" models of Arlandini et al. (1999) and 
Travaglio et al. (1999). 

Key words. Line: formation - Nuclear reactions, nucleosynthesis, abundances - Stars: abundances - Stars: late-type - Galaxy: 
evolution 

1. Introduction elements whose solar system isotopic abundances are domi- 

nated by the ^-process and those due mainly to the r-process. 



The elements heavier than the iron peak are mainly produced ^he europium-to-barium abundance ratio is particularly sen 
^ through neutron capture reactions in two main processes, the ^^^^^^ ^^ whether nucleosynthesis of the heavy elements oc- 
P. , 5-process (slow) and r-process (rapid). As supporte d by m any ^.^^^^ ^y the s- or r-process. Other important indicators are 
observational and theoretical results (Travaglio et aLEH^l and t^e [Eu/La] and [Nd/Ba] abundance ratios. The most recent 
references therein), .v-nuclei are mainly synthesized during the determinations of the [Eu/Ba] abundance ratios in the sam- 
thermally pulsing asymptotic giant branch phase of low-mass pjg^ ^f metal-poor stars were made by McWilUam ( fT998t . 
stars (2 - 4 Mo). The r-process is associated with explosive g^^-is et al. (2000), Fulbright (2000 1, Mashonkina & Gehren 
conditions in SNell (see Thielemann et al.lUnTfor a general ,2000} l2003l . Mishenina & Kovtyukh (l200!t and Bensby et 
review). Following the pioneering observational study of Spite ^j ,2005\. The [Eu/La] abundance ratios in 159 giant and 
& Spite 41978), Truran 11981JI has put forward the suggestion ^warf stars have been measured by Simmerer et al. (2003; 
that most (if not all) neutron capture elements in very metal- j^e [Nd/Ba] abundance ratios in 46 stars have been derived 
poor stars are of r-process origin. Much observational effort ^y Mashonkina et al. (i2004< . Stellar element abundance ra- 
was invested in testing this idea (for review, see Truran et al. ^^^^ ^ave to be compared with theoretical predictions. However, 
12002). To understand the relative importance of the two neu- ^^ present, the nuclesynthesis theory is not able to predict 
tron capture synthesis mechanisms throughout the Galaxy his- j^e yields of the r-process. Based on the hypothesis that the 
tory, stellar element abundance comparisons are made between ^-contributions in the solar system are of primary origin, the 
r-process abundance fractions are obtained by subtracting the 



Send of/print requests to: L. Mashonkina ^-contributions at the epoch of solar formation, estimated with 
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culate the ^-process abundance fractions. The "classical" ap- 
proach relies upon reproducing the product of n-capture cross- 
section and .s-process abundance (the "crN" curve). "Stellar" 
s-process models are based on nucleosynthesis computations 
in low- and intermediate-mass AGB stars. The recent "classi- 
cal" (Arlandini et al. 1999) and "stellar" (Arlandini et al. 1999 
Travaglio et al. .l999j models give consistent results within the 
respective uncertainties for the s-process contributions to the 
most abundant Ba and Eu isotopes in the solar system mat- 
ter, with the only exception being '^^Ba. Its r-residual equals 
14% to 16% in the "stellar" models and 0% in the "classi- 
cal" model. Since '"'^^Ba is the most abundant solar Ba isotope, 
the r-residual of total Ba is larger in the "stellar" model by a 
factor of 2.4 than in the "classical" one. As a result, the "stel- 
lar" and "classical" models lead to significantly different val- 
ues of the solar abundance ratio of Eu to Ba contributed by 
the r-process: taken relative to the total abundances, [Eu/Ba]r 
- 0.69 and [Eu/Ba], = 1.06, respectively. The "stellar" model 
predicts the isotope abundance ratio '^^^Ba : '^^Ba : '-"^Ba = 26 
: 20 : 54 in a pure r-process nucleosynthesis, while no Ba even 
isotope is predicted in the "classical" model. 

To reconstruct the evolutionary history of neutron-rich el- 
ements in the Galaxy it is, thus, very important to disentan- 
gle the fractions of Ba even and odd isotopes and inspect their 
abundance as a function of Galactic age. Isotopic shifts of the 
barium lines are negligible compared to the line width. A de- 
termination of the Ba even-to-odd isotope abundance ratio in 
stars becomes possible due to the significant hyperfine struc- 
ture (HFS) affecting the Ban resonance lines of the odd iso- 
topes. In total, the resonance line AA55A has 15 components 
spread over 58 mA. The larger the fraction of the odd isotopes, 
the stronger the HFS broadening of /14554 and the larger the 
energy absorbed in this line. The total width of the patterns of 
the Ban subordinate lines, /15853, /16141, and /16496, is much 
smaller, 8 mA, 8 mA, and 23 A, correspondingly. 



Cowley & Frey J1989> pointed out the importance of ac- 
counting for HFS in Ba abundance determinations; however, 
the first attempt to estimate the fraction of the odd isotopes 
of Ba in four metal-poor stars was made by Magain & Zhao 
J1993a> . They used the Ban subordinate lines, /15853 and 
/16496, to derive the total Ba abundance. The Ba abundance 
was then deduced from the /14554 line for various isotopic mix- 
tures, the proportion of the odd isotopes being changed until 
agreement with the subordinate lines was obtained. Magain & 
Zhao found an enhancement of the odd isotopes in these stars 
as compared to solar system matter, in agreement with the ex- 
pectations. Using the same approach, Mashonkina et al. ('1999") 
did not find significant deviation of the Ba isotopic mixture in 
the halo star G246-38 ([Fe/H] = -2.20) from the solar one. In 
their recent study, Mashonkina et al. ("2003) show a distinction 
between the different Galactic stellar populations with respect 
to the Ba even-to-odd isotope ratio. The halo stars reveal, on 
average, equal amounts of the odd and even isotopes of Ba; a 
mean ratio 65 : 35 (+10%) was obtained for the sample of thick 
disk stars, whereas the solar ratio 82 : 18 (Anders & Grevesse 
11989 ^ adjusts to observations of the Ban lines in the sample of 
thin disk stars. 



Magain & Zhao ( 11993bl l and Magain (11995') have suggested 
a different method based on measuring the broadening of the 
Ba II A 4554 line. They have applied it to the extremely metal- 
poor star HD 140283 and found the low fraction of the odd iso- 
topes, 8%, that is consistent with a pure ^-process production 
of barium. Using the same method and new high resolution (R 
- 200000), high signal-to-noise ratio (S/N ^ 550) spectra of 
HD 140283, Lambert & Allende Prieto (2002i draw the oppo- 
site conclusion. They find that the r-process mixture of the 
Ba isotopes predicted by the "stellar" model of Arlandini et al. 
([1999) provides a fair fit to the observed Ban A 4554 profile. 

This study is intended to inspect a history of Ba isotopic 
fractions in the Galaxy and, thus, a relative importance of the 
r- and s-process in heavy element production through the 
Galaxy evolution. We determine the abundance ratios of even 
to odd isotopes of Ba in the 25 selected stars representing the 
older (thick disk) and younger (thin disk) Galactic stellar pop- 
ulations. We use the method applied in our previous papers. 
The results are based on the high resolution (~ 60000) high 
signal-to-noise ratio (S/N > 200) spectra and non-local ther- 
modynamical equilibrium (NLTE) line formation. 

The paper is organized as follows. In Sect. |2lwe discuss 
an accuracy of the used atomic parameters, oscillator strengths 
and van der Waals damping constants C6, and check the inves- 
tigated Ban lines with respect to blending lines. We test our 
method on solar Ban lines. Stellar sample, observations, and 
stellar parameters are described in Sect.|3] In Sect.l^we present 
the obtained Ba even-to-odd isotope ratio in stars and investi- 
gate their errors caused by uncertainties of atomic data and stel- 
lar parameters. In Sect.|5l stellar fractions of the odd isotopes 
of Ba are inspected as indicators of the r/i-process nucleosyn- 
thesis, and conclusions are given. At the end, the outlook for a 
determination of Ba isotopic fractions in very metal-poor stars 
is estimated. 



2. Method of calculations 

Homogeneous blanketed model atmospheres computed with 
the MAFAGS code (Fuhrmann et al. 1997^ are used in the anal- 
ysis of both solar and stellar spectra. 

2. 1 . Atomic parameters 

In the metallicity range of our stellar sample, not only the res- 
onance line, but also one or both of the used subordinate lines 
either lie on the damping part of the curve of growth or are 
saturated. In the first turn, our analysis requires very accurate 
atomic parameters. 

As in our previous analysis, in this study we use the 



absolute oscillator strengths from Reader et al. ( I1980I I: 
logg/(^4554) = 0.162, loggf(A5853) = -1.0, and 
log gf {A6496) - -0.377. The most recent measurements 
of Davidson et al. jTW^ give log g/(^ 4554) = 0.140, 
log g/(^ 5853) = -0.91, and log gf (A6496) = -0.407. We 
show below that applying the data of Davidson et al., we can- 
not agree when fitting parameters of the solar Ban A 6496 and 
/i5853 lines. We do not use the third line of the Ban 5d - 6p 
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multiple t, /1 6141 which is blended with the strong Fel line 

(Fig.GJ. 

In our previous studies, we used the data on hyperfine 
structure components of the Ba ii A 4554 calculated by Biehl 



J1976> on the base of Brix & Kopfermann (1952i measure- 
ments. Having reviewed the literature for more recent data, we 
found that the new fine structure constants of '^^Ba and '^^Ba 
obtained by Blatt & Werth ( 1982 1 and Becker & Werth ( 1983) 
lead to a negligible change of the wavelength separations by 
no more than 1 % for the HFS components of the resonance 
line. In order to have common input data with other authors 
that makes a direct comparison of the results possible, in this 
study all calculations of /1 4554 are performed using the HFS 
pattern published by Mc William ( 1998 1. We emphasize that the 
line profiles computed with the data of Biehl and Mc William 
become identical after their relative shift by 0.001 A. Tabled 
contains the list of the HFS components of the Ban /1 4554. 
The product of fij and fractional isotope abundance e is given 
for the Ba isotope mixture in the solar system matter, '■''^Ba 
: i^^^Ba : '^^Ba : '"Ba : '^'^Ba = 2.4 : 6.6 : 7.9 : 11.2 : 71.8 
(Anders & Grevesse 1989 1 and that for a pure r-process pro- 
ducti on of Ba, '^^Ba : '"Ba : '^'^Ba = 26 : 20 ; 54 (Arlandini 
et al. n^y^ . Total fij = 0.727. The 23 mA wide patterns of the 
/16496 line are considered as a four component model accord- 
ing to Rutten ( 1978 1. Variation in fractional isotope abundances 
between the solar and pure r-process cases has negligible ef- 
fect on the /16496 line profile, and we use /xe computed for the 
solar mixture of Ba isotopes. We neglect the hyperfine structure 
of/15853. 

The investigated Ban lines are strongly affected by van der 
Waals damping. We apply the van der Waals damping con- 
stant logCe = -31.65 to /1 4554. If logCe is fixed, solar Ba 
abundance is immediately determined from the fitting of the 
solar /14554 line wings. Using the MAFAGS solar model at- 
mosphere, we find logfioBa = 2.21 '. This value agrees within 
error bars with the meteoritic Ba abundance given by Anders 
& Grevesse ( 119891 logemet,Ba - 2.21 + 0.03) and Grevesse 

2.22 + 0.02). Recent determina- 



et al. J 19961 logemet.Ba 
tions by Asplund et al. J2005t lead to the smaller meteoritic 
(log£n,et,Ba = 2.16 + 0.03) and solar (logSoBa = 2.17 + 0.07) 
Ba abundance. For comments, see Note added in proofs. 

The Van der Waals damping constant log C(, of the subordi- 
nate lines was obtained empirically from the fitting of the solar 
line profiles under the requirement that unique set of fitting pa- 
rameters (Ba abundance, microturbulence value, and C^ value) 
must adjust observations of the spectral lines belonging to the 
same multiplet, 5d - 6p. Far wings of /1 6496 impose the up- 
per limit of log Cfi, provided that solar Ba abundance is given. 
We adopt logSoBa - 2.21. In this case, logCe of /1 6496 can- 
not be larger than -31.3. Applying this value to /1 5853 and 
based on NLTE analysis we obtain a microturbulence value 
Vmic - 0.85 kms"'. The best fit of each line of the multiplet 
5d - 6p is achieved at log SBa = 2.21, log Ce = -31.3, and V^ic 
= 0.85 kms"' (Fig-[D- In Sect. 0] we will discuss an effect of 
variation in log Ce on Ba isotope ratio in stars. Analysis of the 
solar /1 6496 and /1 5853 line profiles was also performed with 



Table 1. Atomic data for the HFS components of the Ban 
/1 4554. / is the relative strength, the product of oscillator 
strength and fractional isotope abundance, /xe, is given for two 
Ba isotope mixtures corresponding to the solar system matter 
(column "solar") and r-process yields predicted by Arlandini 
etal. (1999) 



A, A 


/ 


/X6 




solar 


r-process 


"^Ba 








4554.032 


1.000 


0.0176 


- 


'-"Ba 








4554.003 


0.1562 


0.0075 


0.0292 


4554.004 


0.1562 


0.0075 


0.0292 


4554.005 


0.0625 


0.0030 


0.0117 


4554.051 


0.4375 


0.0209 


0.0818 


4554.054 


0.1562 


0.0075 


0.0292 


4554.055 


0.0313 


0.0015 


0.0058 


'36Ba 








4554.033 


1.000 


0.0568 


- 


'"Ba 








4553.999 


0.1562 


0.0129 


0.0232 


4554.001 


0.1562 


0.0129 


0.0232 


4554.002 


0.0625 


0.0052 


0.0093 


4554.054 


0.4375 


0.0360 


0.0650 


4554.056 


0.1562 


0.0129 


0.0232 


4554.057 


0.0313 


0.0025 


0.0046 


'38Ba 








4554.034 


1.000 


0.521 


0.392 



We refer to abundances on the usual scale where log en = 12. 



the oscillator strengths from Davidson et al. ( I1992> . Solar Ba 
abundance was fixed at logSoBa = 2.21. Fitting parameters of 
two lines turn out to be different: A 6496 gives logCe = -31.2 
and Vmic = 0.9 kms ', while the significantly lower microtur- 
bulence value Vmic = 0.7 kms"' is required to fit the /1 5853 
line profile using log €(, = -31.2. 

The Ban A 4554.034 A Hne is blended by the Cr I 
.14553.945 and Zrn ^4553.934 Unes in the blue Hne wing 
and by the CH molecular line at 4554.142 A in the red 
line wing (Fig. O. We have requested the Vienna Atomic 
Line Data base (VALD, Kupka et al. l^W) and the NIST 
atomic spectra database (http://physics.nist.gov/PhysRefData) 
to search for atomic parameters of blending atomic lines. The 
same value \og gf = -0.73 based on measurements of Wujec 
& Weniger (1981) is given in NIST and VALD for the Cr I 
.14553.945. NIST does not contain the Zrii .14553.934 Hne, 
but VALD provides \og gf - -0.57, taken from Cowley & 
Corliss ( 1983l. for this line. The oscillator strength of the 
CH /l 4554. 142 molecular line was fitted to reproduce the ob- 
served blend profile; we have obtained \og gf - -3.2. The 
Ban .15853.675 A line is blended by the Fe 1 .1 5853. 682 Hne 
(Fig.Hli. Based on calculations of Kurucz & Bell ( IT9951 . VALD 
gives log gf (Fe I .15853.682) = -2.148. Recently, Kurucz 
(http://cfaku5.cfa.harvard.edui calculated the new value \og gf 
- -231 \ for this line. We note that the predicted oscillator 
strength of Fe I /1 5853 .682 has been reduced by a factor of 2.3, 
compared to the first estimate of Kurucz & Peytremann \\915\ . 
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log gf — -2.0. This line is absent in the NIST database. The 
Ban /16141.713 A is strongly blended by the Fe I ^6141.732 
line and not used in further analysis of stellar spectra. With 
\oggf (Fe 1/16141) - -1.61 taken from the NIST database, we 
can achieve a good fitting of the solar blend profile (Fig. ^. 
VALD gives the larger oscillator strength, \og gf - -1.459, for 
Fe I /16141. The contribution of blending lines is treated us- 
ing the LTE assumption. The chemical abundances are taken 
from Anders & Grevesse ( 1989 1 in the solar case, and we as- 
sume that in the investigated stars, Cr abundance follows the 
iron one with [Cr/Fe] - 0, and Zr abundance follows the bar- 
ium one with [Zr/Ba] = 0. 




2.2. NLTE line formation 

Our results are based on NLTE line formation for Ban. The 
method of NLTE calculations was developed and described ear- 
lier (Mashonkina & Bikmaev [19961 Mashonkina et al. 1999). 
We use a revised version of the DETAIL program (Butler 
& Giddings 1985) based on the accelerated lambda iteration 
method to solve the coupled radiative transfer and statisti- 
cal equilibrium equations and the SIU code (www.usm.uni- 
muenchen.de/people/reetz/siu.html) to compute the synthetic 
line profiles. 

In Fig.n we present the best fits of the observed Ban sub- 
ordinate line profiles in the solar flux atlas of Kurucz et al. 
<1984> achieved using the NLTE approach and LTE assump- 
tion. We emphasize that the unique value of Ba abundance, 
log eBa - 2.21, adjusts observations of all these lines. The same 
value, logCe - -31.3, is used for each line in both the NLTE 
and LTE calculations. The same microturbulence value, V^^ic 
= 0.85 kms \ is required when the NLTE approach is ap- 
plied while different values in the LTE case are used: Vmic = 
0.95 kms"^ for .i5853, V,™c = 1-05 kms"' for .16496, and 
Vmic = 1 kms"^ for /16141. Our synthetic flux profiles are 
convolved with a profile that combines a rotational broaden- 
ing of 1.8 kms"' and broadening by macroturbulence with 
a radial-tangential profile. The macroturbulence values Vmac 
= 3.1 kms"' and 3.0 kms"' are required for /1 5853 in the 
NLTE and LTE cases, respectively. The corresponding values 
for /1 6496 are 2.3 kms ' and 2.1 kms ' and for /i6141, in 
both cases, Vmac - 2.9 km s"' . 



2.3. Can we derive tine solar Ba isotope abundance 
ratio? 

Assuming the solar system matter Ba isotopic fractions '•'"^Ba : 
'^^Ba : '^"^Ba : '^^Ba : '^'^Ba = 2.4 : 6.6 : 7.9 : 11.2 : 71.8, total 
Ba abundance log SBa - 2.21, log Ce = -31.65, and using Vmic 
= 1 kms"', we can achieve a good fitting of the solar /1 4554 
profile. In Fig. |2l we show the best fit obtained by applying 
Vmac = 3.5 kms"'. 

The value Vn-iK - 1 km s"' obtained from the resonance line 
does not agree with the value Vmic - 0.85 kms"' found from 
the subordinate lines. We are aware that a phenomenological 
description of the velocity field using a microturbulence value 
is a rough approximation and that it fails to treat spectral lines 




6496.70 6496.80 6496.90 6497.00 6497.10 




6141.50 6141.60 6141.70 6141.80 6141.90 6142.00 



Fig. 1. The best fits of the solar flux Ba II subordinate line pro- 
files (Kurucz et al. 1984, bold dots) achieved using the NLTE 
approach (continuous line) and LTE assumption (dashed line). 
Solar Ba abundance log SBa = 2.21; log Ce = -31.3. NLTE pro- 
files: each is calculated using Vmic - 0.85 kms"'. In the LTE 
case, y^ic = 0.95 kms"' (.15853), 1.05 kms"' (.16496), and 
1 kms"' (/16141). Blending lines are shown in each panel by 
a dash-dotted line. See text for more details. 



formed over a wide range of depths. The A 4554 line core is 
formed in very outer layers near log T5000 = -4.7, while radia- 
tion of the subordinate lines comes from the deeper layers, be- 
low log T5000 = -2.8 for A 6496 and below log rsooo = -1.9 for 
A 5853. We have tried to answer the question of whether we are 
able to derive a solar atmosphere Ba isotope abundance ratio 
from the HFS affecting the A 4554 line profile without a refer- 
ence to the subordinate lines. In addition to the even-to-odd iso- 
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4553.70 4553.80 4553.90 4554.00 4554.10 4554.20 4554.30 



Fig. 2. Synthetic NLTE flux profile of the Ba II resonance line 
corresponding to the Ba even-to-odd isotope ratio of 82 : 18 
(continuous line) compared with the solar flux profile (bold 
dots, Kurucz et al. 1984). The diff'erences between observed 
and calculated spectra, (O - C), multiplied by a factor of 5, are 
shown for various even-to-odd isotope ratios in the lower part 
of figure. Blending lines are shown by a dash-dotted line. See 
text for more details 



tope ratio 82 : 18 corresponding to solar system matter, a fitting 
of the solar A 4554 line profile was made assuming the ratios 89 
: 11 and 75 : 25. We suppose that the proportion '■'^Ba : '^^Ba 
does not change. It is not the case if the r/ i-process contribu- 
tion to Ba varies. However, as can be seen from Tabled the 
/1 4554 line profile depends on the total fractional abundance 
of the odd isotopes rather than on the proportion '''''Ba : '''^Ba. 
Our test computations for the isotope mixture '^^Ba : '^^Ba : 
(i34gjj + i36g^ + i38g^-) = 2 : 9 : 89 predicted by Arlandini 
et al. J1999> for a pure .s-process production of Ba give the 
/14554 hne profile fully consistent with that for the ratio '^^Ba 
: ™Ba : C^^Ba + '^^Ba + '^'^Ba) = 5.5 : 5.5 : 89. In calculations 
with various even-to-odd isotope ratio, a microturbulence value 
was a free parameter, and in each case the best fit was found 
from the minimum root mean square (rms) difiference between 
observed and calculated spectra. Vmac was fixed at 3.5 kms"'. 
The minimum rms is achieved at Vmic = 1.15 kms"' for the 
even-to-odd isotope ratio 89 : 11, and Vmic = 0.9 kms"' for 
75 : 25. In the latter case, a shift of the observed profile by 
-0.0008 A is required. Fig. |3 shows the obtained differences 
between observed and calculated spectra, (O - C), multiplied 
by a factor of 5, for various isotope ratios. 

The rms values quoted in Fig. |2l are very similar for vari- 
ous isotope mixtures. The obtained Vmic values 0.9 kms"' to 
1.15 kms"' are within the interval of widely used solar mi- 
croturbulence values, and no value can be preferred. The (O - 
C) values are the largest in the line core within +0.07 A from 
the line center, independent of isotope ratio. The A 4554 line 
core is most probably influenced by a non-thermal and depth- 
dependent chromospheric velocity field that is not part of our 
solar model. We note that the chromospheric temperature rise, 
as well as the horizontal temperature inhomogeneity seen as 
solar granulation, is expected to be less important because Ban 
is the dominant ionization stage. 



We conclude from these results that a precise determina- 
tion of the proportion of the even and odd isotopes of Ba in 
the solar atmosphere matter requires a correct treatment of the 
atmospheric velocity field that is only possible on the basis of 
hydrodynamic calculations. In metal-poor stars, the A4554 line 
forms in the deeper atmospheric layers compared to the Sun 
and in smaller extension is expected to be aff'ected by hydrody- 
namic phenomena of the very surface layers. 

3. Stellar sample, observations, and stellar 
parameters 

All the 25 stars were selected from Fuhrmann's f2004^ list. 
Each star was observed with a resolution of ~ 60 000, at least 
twice, and for each star both the Ban subordinate lines, /15853 
and /16496, are available in its spectra. Spectra were obtained 
by Klaus Fuhrmann using the fiber optics Cassegrain echelle 
spectrograph FOCES at the 2.2m telescope of the Calar Alto 
Observatory in 1997 - 2000. The signal-to-noise ratio is 200 or 
higher in the spectral range A > 4500 A. Typical line profiles 
are seen in Figs.|3and|5] 

The list of stars under investigation is given in Table |2 
Stellar parameters have been determined spectroscopically by 
Fuhrmann (2004i. Their errors are estimated as Aretf - 80 K, 
A log ^ = 0.1, A[Fe/H] = 0.1 dex, and Ay^c = 0.2 kms"'. We 
use also his identification of a membership of individual stars 
in the particular stellar population that is based on the star's 
kinematics, metallicity, a-enhancement, and age. Our sample 
includes 15 thick disk stars, 9 thin disk stars, and one halo star. 

4. Results 

4.1. Stellar total Ba abundances 

Ba abundances were derived from the Ban /1 5853 and /1 6496 
line profile fitting, and then the average value was calculated. 
The [Ba/H] values calculated using solar Ba abundance logSBa 
- 2.21 are presented in Table |3 We note that in this paper we 
revise our earlier determinations (Mashonkina & Gehren l2001> 
of Ba abundance in the investigated stars. Contrary to the pre- 
vious analysis, we use the same value log Cg for both subordi- 
nate lines and treat the /1 5853 blend taking into account both 
the Ban and Fe I lines (see Sect. l2.lt . A reduction of the van 
der Waals damping constant from logCe = -30.6 for /1 5853 
and log Ca = -3 1 .2 for A 6496 used in 2001 to our present esti- 
mate log Ce = -3 1 .3 (Sect. 12. H leads to an increase of derived 
Ba abundance, while accounting for absorption by the blending 
of the Fe I line leads to a decrease of Ba abundance from the 
/i5853 blend. As a result, in the two most metal-poor stars of 
our sample (HD 103095 and HD 64606), the revised Ba abun- 
dance is smaller (van der Waals broadening of the subordinate 
lines is weak and blending of Ban /i 5853 prevails) by 0.01 dex 
compared to our previous data. Also in the remaining stars, the 
revised Ba abundance has become larger by 0.02 - 0.04 dex, 
with the only exception being that HD 1 17176 now reveals Ba 
abundance that is higher by 0.06 dex. The latter star shows 
the largest discrepancy in logSBa derived from the /1 5853 and 
/1 6496 lines. 



L. Mashonkina and G. Zhao: Barium even-to-odd isotope abundance ratios in thick disk and thin disk stars 



Table 2. Stellar parameters, Ba abundances, and the fractions 
of the odd isotopes of Ba (column odd, in %) of the selected 
sample. The column low - up contains the lower and upper lim- 
its of that fraction. In the column "Note", the notations 0, 1, and 
2 refer to the thin disk, thick disk, and halo stars, respectively. 
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Ba abundance determined from the /1 5853 line turns out 
to be systematically lower than that from the second line; the 
Ba abundance is lower by -0.03 + 0.02 dex in the stars with 
[Ba/H] > -0.35, by -0.01 ±0.03 dex in the more Ba-poor stars, 
and by -0.02 + 0.03 dex for the whole sample (Fig. 0] top 
panel). The abundance difference loge(/l5853) - loge(/i6496) 
shows no correlation with T^.^ (Fig. 0] bottom panel). Ba 
abundance from the /1 5853 line can be underestimated due 
to overestimation of the oscillator strength of the blending 
Fe I /15853 line. Using \og gf - -2.371 calculated recently 
by Kurucz I|httpy/cfaku5 .cfa.harvard.edut, we derive the larger 
Ba abundance from the /1 5853 blend: it is larger by 0.01 dex 
in HD9407 ([Ba/H] = -0.01) and by 0.02 dex in HD 68017 
([Ba/H] = -0.50) and HD 64606 ([Ba/H] = -1.00). Another 
source of the obtained discrepancy between the two Ban lines 
is the uncertainty of the van der Waals damping constant of 



0.05 - 0.18 dex according to Barklem J2006t . In the next sub- 
section, we will discuss an effect of total Ba abundance error 
on the Ba even-to-odd isotope ratio in stars. 
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Fig. 3. Synthetic NLTE (continuous line) flux profiles of the 
Ba II subordinate lines compared with the observed FOCES 
spectra (bold dots) of HD 10519. The same value [Ba/H] - - 
0.69 is obtained from both lines 
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Fig. 4. NLTE abundance differences loge(/l5853) - 
log e(/l6496) versus total Ba abundance [Ba/H] (top panel) and 
Teff (bottom panel) 

4.2. Ba even-to-odd isotope ratio in stars 

The fraction of the odd isotopes of Ba is derived from the re- 
quirement that Ba abundances derived from the Ban resonance 
line and the subordinate lines must be equal. We emphasize 
that absolute Ba abundance log SBa is used in this procedure, 
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Fig. 5. Synthetic NLTE flux profile corresponding to log SBa = 
1.52 and Ba even-to-odd isotope ratio 58 : 42 (continuous line) 
compared to the observed profile (bold dots) of the Ba II res- 
onance line in HD 10519. logSBa - 1-51 for the isotope ratio 
54 : 46, logSBa = 1-54 for 65 : 35, and logeea = 1-64 for 82 
: 18. The (O - C) values multiplied by a factor of 5 are shown 
for various isotope mixtures in the lower part of the figure. See 
text for more details 



and, thus, the obtained fraction does not depend on the adopted 
value of solar Ba abundance. 

Figure [S] illustrates a determination procedure for one of 
the stars of our sample, HD 10519. The same Ba abundance, 
logfiBa = 1.52, is obtained from both subordinate lines. The 
A 4554 line gives log SBa - 1 -64 at the assumption of the iso- 
tope mixture in solar system matter and log eBa =1.51 for the 
pure r-process yields predicted by Arlandini et al. ( I1999t . As 
noted above, the synthetic A 4554 line profile depends on the to- 
tal fractional abundance of the odd isotopes rather than on the 
exact abundance of each odd isotope. Therefore, we vary the 
even-to-odd isotope ratio (from here on, "the isotope ratio") in 
calculations of /1 4554 until the obtained Ba abundance agrees 
with that from the subordinate lines. The required Ba abun- 
dance logfiBa - 1-52 is achieved from /1 4554 when the isotope 
ratio equals 58 : 42. The following points should be noted. 

- When we go from the isotope ratio 82 : 18 to 54 : 46, the 
relative shift of the observed profile by 0.001 A is required 
to achieve the best fit. 

- We define the best fit from a requirement of the minimum 
rms difference between observed and calculated spectra. 
The equivalent width corresponding to the best fit may be 
different for different isotope ratios, e.g., in our calcula- 
tions for HD 10519, the equivalent width measured over 
the same spectral interval, 4553. 8A- 4554.3A, W(/14554) 
= 154.1 mA for the isotope ratio 82 : 18 and 152.2 mA for 
54 : 46. This can be understood. The smaller fraction of the 
odd isotopes that the larger Ba abundance is required to fit, 
in the observed line profile, and the calculated line wings, 
which never perfectly fit the observed ones (Fig. |5|l, con- 
tribute to W,) in the greater extension, when compared to 
the case of the larger fraction of the odd isotopes with the 
smaller Ba abundance. 



- The macroturbulence value was allowed to be free because 
we have no arguments to fix it. An increase of the frac- 
tion of the odd isotopes cancels, in part, a saturation of 
/14554, and the line formation depth shifts downward; e.g., 
in our calculations for HD 10519, a macroturbulence value 
grows from 4. 25km s"' to 4.5km s"' when the isotope ratio 
changes from 82 : 18 to 54 : 46. 

Uncertainty of the fractional abundance of the odd isotopes 
ofBa in stars. Random errors of a desirable fraction are due to 
total Ba abundance error and the uncertainty of stellar parame- 
ters, Teff, logg, and microturbulence value. Systematical errors 
are caused by the uncertainty of the Ba ii A 4554 atomic param- 
eters, namely, fj and Cf, values. Varying all mentioned param- 
eters, we have performed the test calculations for the three rep- 
resentative stars of our sample with different values of [Ba/H], 
HD9407 ([Ba/H] = -0.01), HD 68017 ([Ba/H] = -0.50), and 
HD 64606 ([Ba/H] = -1.00). The effect on the Ba odd isotope 
fraction in these stars is shown in Table|3l 

The most important source of total Ba abundance error in 
the stars with metallicity close to solar one is the uncertainty of 
the van der Waals damping constant of the subordinate lines. 
A variation in log Ce of 0. 1 translates to the 0.03 dex variation 
in logSBa- When /1 4554 is saturated, a reduction of Ba abun- 
dance obtained from the subordinate lines by 0.01 dex requires 
an increase of the fractional abundance of the odd isotopes by 
~2% to ~3% depending on the resonance line strength. A vari- 
ation of 80 K in Teff translates to the uncertainty of the obtained 
fraction of 3% to 4%. In the stars with [Ba/H] < -0.5, a vari- 
ation in Vmic has nearly the same effect on both the resonance 
and subordinate lines, and no change of the Ba odd isotope 
fraction is required. In the stars with metallicity close to solar 
one, the resonance line lies on the damping part of the curve of 
growth and is weakly sensitive to y^ic variations, while this is 
not the case for the subordinate lines. The uncertainties of total 
Ba abundance, T^s, logg, and Vmic, in total, result in a random 
error of the obtained fraction of 7%, 9%, and 10% in the stars 
with [Ba/H] - 0, -0.5, and -1, correspondingly. There can be 
also a systematical difference of 4% to 5% between the thick 
disk and thin disk stars, due to a discrepancy in the mean value 
log eBa('^5853) - log eBa(^6496) between the samples of thick 
disk and thin disk stars. 



Table 3. Effect on the Ba odd isotope fraction (in %) caused by 
uncertainties of atomic data and stellar parameters 
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A variation of 0. 1 dex in the van der Waals damping con- 
stant of the A 4554 line produces a strong effect, 7 % to 8 % 
of the derived fraction of the odd isotopes of Ba; however, this 
is a systematical effect. It depends only weakly on Ba/H and 
therefore will not affect a difference in the fraction of the odd 
isotopes between different stars. 

We suppose that in each individual star, the spread of Ba 
abundance derived from two subordinate lines fully reflects 
both random and systematical errors of total Ba abundance, in- 
cluding effects of the uncertainties of atomic parameters (gf 
and C(, values), microturbulence value, and NLTE treatment. It 
is worth while to remember that due to different strengths of 
the Ban ^6496 and ^5853 lines (/(/l6496)//(/l5853) ^ 4), 
they lie in each star on different parts of the curve of growth. 
The uncertainty of a desirable fraction in each star is evaluated, 
taking into account its individual Ba abundance error and un- 
certainty of the /14554 line profile fitting. In terms of Ba abun- 
dance, the latter value is estimated to be 0.01 dex, based on the 
analysis of the (O - C) values. For example, in HD 10519, both 
subordinate lines give the same log SBa and the lower and up- 
per limits of the odd isotope fraction are calculated taking into 
account only the uncertainty of the A 4554 line profile fitting 
(0.01 dex). In the range of the odd isotope fraction between 
35% and 46%, characteristic of this star, a Ba abundance vari- 
ation of 0.01 dex produces the 4 % variation in the fractional 
abundance of the odd isotopes. Therefore, our estimate of the 
fraction of the odd isotopes of Ba in this star is 42 ± 4%. 

The obtained fraction of the odd isotopes, with its lower 
and upper limits in the stars of our sample, are presented in 
Tableia 

5. Stellar Ba odd isotopes as an indicator of the 
r/i^-process nucleosynthesis 

In Fig. |6l we plot the fractional abundance of the odd iso- 
topes of Ba versus the total Ba abundance and the [Eu/Ba] and 
[Eu/Fe] abundance ratios. Based on the spread of Ba abundance 
in each star (Fig.|3}, error bars of the [Ba/H] value do not ex- 
ceed 0. 1 dex. The [Eu/Fe] abundance ratios are taken from our 
earlier determinations (Mashonkina & Gehren 1200 lb . where 
their errors are estimated as A[Eu/Fe] = 0.05 dex. We note 
that both Eu and Fe abundances are derived from weak spectral 
lines of the dominant ionization stage, Eu ii and Fe ii, respec- 
tively, and the uncertainty of stellar parameters has a negligible 
effect on the derived Eu/Fe abundance ratio. Error bars of the 
[Eu/Ba] value are mainly defined by errors of Ba abundance. 

It is clearly seen that the obtained fraction grows towards 
the lower Ba abundance. Its value in the thick disk stars ranges 
between 27% and 42%, with the mean value 33+4%. This indi- 
cates the higher contribution of the r-process to barium in the 
thick disk stars compared to the solar system matter in agree- 
ment with our expectations. Analysis of the [ff/Fe] abundance 
ratios (we cite here only the first such studies: Gratton et al. 
[19961 EUnni FuhrmannTWT Prochaska et al.lUnO") suggests 
that the thick disk stellar population is nearly the same age as 
the Galactic halo one. This conclusion was confirmed from an 
analysis of Eu/Ba and Nd/Ba abundance ratios (Mashonkina & 
Gehren l2000t Mashonkina et al. 2004. Bensby et al. 2005.). The 
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Fig. 6. The fractional abundance of the odd isotopes of Ba ver- 
sus [Ba/H] (top panel), [Eu/Ba] (middle panel), and [Eu/Fe] 
(bottom panel). The uncertainty of individual values is shown 
by vertical lines. Horizontal lines indicate the solar fraction of 
the odd isotopes, 18%, and the values predicted by Arlandini 
et al. (1999) and Travaglio et al. (1999) for a pure r-process, 
46%, and pure i-process, 11%, production of barium 



significant overabundance of Eu and Nd relative to Ba found 
in the thick disk stars at [Fe/H] < -0.3 says that these stars 
formed in the early Galaxy when SNell dominated the synthe- 
sis of heavy elements. 

Both the [Eu/Ba] abundance ratio and the fractional abun- 
dance of the odd isotopes of Ba are sensitive to the relative con- 
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tribution of the r-process to heavy element synthesis, and it is 
natural to expect their correlation. Figure|6lshows that such cor- 
relation really exists. It is more clearly seen in the plane frac- 
tional abundance of the odd isotopes of Ba - [Eu/Fe] because 
the range of [Eu/Fe] values in the thick disk stars is larger than 
that of [Eu/Ba], and errors of the [Eu/Fe] values are smaller 
thanthatof [Eu/Ba]. 

The obtained data provide observational constraints to 
models of the s- and r-process nucleosynthesis. The "classi- 
cal" model of the i-process based on the newest accurate mea- 
surements of the stellar n-capture cross-sections (Arlandini et 
al.[T999) predicts a pure i-process origin of solar even isotopes 
of Ba. In this model, no even isotope of Ba appears in stars 
formed in the early Galaxy before the onset of the i-process 
in AGB stars. In contrast to this prediction, we find a signifi- 
cant fraction of the even isotopes of Ba, ~ 67%, in old Galactic 
stars, the thick disk stars. The "stellar" models of Arlandini 
et al. 11999), based on stellar AGB models of 1.5 Mq and 
3 Mq with half solar metal abundance, and Travaglio et al. 
<1999> . who integrate s-abundances from different generations 
of AGB stars, i.e., consider the whole range of Galactic metal- 
licities, both predict that solar abundance of isotope '^**Ba is 
contributed to by the s-process only in part: 86% and 84% , re- 
spectively. As a result, in the oldest stars of the Galaxy, formed 
at the epoch of SNell dominance in nucleosynthesis, a frac- 
tion of the even isotopes of Ba is expected to be at the level of 
54%. Ba isotopic fractions found in the thick disk stars favor 
the "stellar" model of heavy element synthesis. 

Acknowledgements. We are very grateful to Klaus Fuhrmann for 
providing reduced stellar spectra and Thomas Gehren for provid- 
ing a Windows version of the code DETAIL. ML acknowledges 
with gratitude the National Astronomical Observatories of Chinese 
Academy of Science for warm hospitality during a productive stay 
in October - December of 2005. We thank the anonymous referee 
for useful remarks and comments. This research was supported by 
the Deutsche Forschungsgemeinschaft with grant 436 RUS 17, the 
Russian Foundation for Basic Research with grant 05-02-39005- 
GFEN-a, the Natural Science Foundation of China with grants NSFC 
10433010 and 10521001, the RF President with a grant on Leading 
Scientific Schools 1789.2003.2, and the Presidium RAS Programme 
"Origin and evolution of stars and the Galaxy". 

Appendix A: The outlook for a determination of Ba 
isotopic fractions in very metal-poor stars 

It would be very important to extend our study to the older and 
more metal-poor stars revealing a pure r-process nucleosyn- 
thesis and to determine from observations the relative yields of 
the even and odd isotopes of Ba in the r-process. At [Ba/H] 
< -2, our method fails to give a reliable value of the fractional 
abundance of the odd isotopes of Ba due to (i) less sensitivity 
of the A 4554 line to its variation, and (ii) less accuracy of total 
Ba abundance. An increase of the fraction of the odd isotopes 
of Ba from 18% to 46% leads to a decrease of Ba abundance 
derived from /i4554 by 0.07 dex at [Ba/H] = -2 and by only 
0.02 dex at [Ba/H] = -3.3. At the same time, the weaker sub- 
ordinate line /1 5853 becomes, in fact, unmeasurable in stellar 
spectra. 
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Fig.A.l. The differences between "observed" and calculated 
spectra, ("O" - C) in "star" 1 (top panel) and "star" 2 (bottom 
panel) for the even-to-odd isotope ratios 82 : 18 and 65 : 35. In 
the latter case, the differences shifted down by 0.05 are shown. 
The HFS components of the odd isotopes appear in two groups 
shifted relative to the resonance line of the even isotopes at 
4554.034 A by 18 mA and -34 mA. Their positions are indi- 
cated by vertical hnes. See text for more details 



Can we determine a desirable fraction based on only the 
HFS broadening of the /1 4554 line without knowing total Ba 
abundance? Let's imagine that we have a perfect observed pro- 
file of A 4554. Then we fit it varying the fractional abundance 
of the odd isotopes of Ba. We have modelled such situation for 
two "stars" using two model atmospheres. 

"Star" I: T^s = 5710 K, \ogg = 4.00, [M/H] = -0.64, 
Vmic = 1.1 kms ', [Ba/H] = -0.69. The theoretical line profile 
of A 4554 calculated assuming the Ba odd isotope fraction of 
42% and convolved with the Gaussian of 4 km s ' and radial- 
tangential profile of Vmac = 4.5 kms"' will serve as the "ob- 
served" line profile. Its W^ = 152 mA. 

"Star" 2: T^n = 6350 K, \ogg = 4.03, [M/H] = -2.07, 
Vmic =1.7 kms"', [Ba/H] = -2.14. The "observed" profile 
of A 4554 was obtained assuming the Ba odd isotope fraction 
of 46% and applying the Gaussian of 2 kms"' and Vmnc - 
4.7 kms"'. Its W^= 50 mA. 

Assuming the fraction of the odd isotope of Ba, odd = 18% 
and then 35%, we obtained the best fits for each case and for 
both "stars". Ba abundance and Vmnc value were allowed to 
vary. The Gaussian was fixed. In practice, the Vmac values turn 
out to be within the uncertainty of determination: Vmac = 4.25 
kms"' when odd = 18% and 4.4 kms"' when odd -35% in 
"star" 1; Vmac = 5.1 kms"' when odd = 18% and 5.0 kms"' 
when odd = 35% in "star" 2. The ("O" - C) differences are 
shown in Fig. lA.ll 

It can be seen from Fig. lA.ll that an accuracy of line profile 
measurements has to be much higher than 0.05% if one wishes 
to determine a desirable fraction within the uncertainty of 10%. 
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Note added in proofs. 

For A 4554, the van der Waals damping constant based on 
theoretical predictions of Barklem & O'Mara (I1998> is equal to 
log C(, =-31 .46 with the uncertainty of 0.05 - 0. 1 8 dex accord- 



ing to Barklem & Aspelund-Johansson (I2005> and Barklem 
([2006 1. Using this value and the MAFAGS solar model at- 
mosphere, we find solar Ba abundance logSoBa = 2.17, in 
well agreement with the recent determinations of the meteoritic 
(logemet,Ba = 2.16 + 0.03) and solar (log SoBa = 2.17 ±0.07) Ba 
abundance by Asplund et al. ( 2005 1. Assuming log C^ =-31 .46 
for A 4554, we checked other stars of our sample with stellar pa- 
rameters close to solar ones, HD 9407, HD 10697, HD 134987, 
and HD 168009 (see Sect.|3, and found systematically lower 
Ba abundance from the Ban resonance line compared to that 
from the subordinate lines, by 0.06 dex to 0.10 dex. Variation 
of stellar parameters and microturbulence value does not help 
to achieve agreement between different lines, and the only way 
out is to reduce collisional broadening of the resonance line. 
Therefore, we apply in this study the value log C(, =-31 .65 to 
A 4554. For the Ba ii multiplet 5d-6p, the van der Waals damp- 
ing constant logCe = -31.3 found empirically in Sect. l2.1l is 
very close to log Ce = -3 1 .28 based on the theory of Barklem 
& O'Mara (1998 1 and accessible via the Vienna Atomic Line 
Data base (VALD, Kupka et al]T^^ 



